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ABSTRACT 

We probe the feasibility of high-frequency radio observations of very rapid flux variations in compact 
active galactic nuclei (AGN). Our study assumes observations at 230 GHz with a small 6- meter class 
observatory, using the SNU Radio Astronomical Observatory (SRAO) as example. We find that 33 
radio-bright sources are observable with signal-to-noise ratios larger than ten. We derive statistical 
detection limits via exhaustive Monte Carlo simulations assuming (a) periodic, and (b) episodic flaring 
flux variations on time-scales as small as tens of minutes. We conclude that a wide range of flux 
variations is observable. This makes high-frequency radio observations - even with small observatories 
- a powerful probe of AGN intra-day variability; especially, those observations complement observations 
at lower radio frequencies with larger observatories like the Korean VLBI Network (KVN). 

Key words : Galaxies: active — Radiation mechanisms: non-thermal — Methods: statistical 



1. INTRODUCTION 

Active galactic nuclei (AGN; for a review, see e.g. 
Beckmann & Shrader 2012) show strong flux varia- 
tions on various timescales. On the shortest timescales, 
intra-day variability (IDV) of radio fluxes was first re- 
ported by Witzel et al. (1986). Since then, such a rapid 
variability has been observed throughout the electro- 
magnetic spectrum mainly in compact objects, such as 
flat spectrum radio quasars (FSRQs) and BL Lacer- 
tae objects (BL Lacs). At frequencies around 5 GHz, 
about half of all FSRQs have been reported to show 
IDV (Lovell et al. 2008). 

A priori, such a rapid variability is hard to under- 
stand. For a given observed variability time scale i var 
the spatial extension of the emission region is limited 
to ct val /(z + 1), with c denoting the speed of light and 
z being the redshift. The observed brightness tem- 
perature is related to the variability time scale like 
Tb oc [ivar [z + 1)] 2 ; observed values are as high as 
T b « 10 18 K (Wagner & Witzel 1995; Fuhrmann et al. 
2008). However, cooling by inverse Compton scattering 
limits the brightness temperature of realistic plasmas 
to T b <; 10 12 K (Kellermann & Pauliny-Toth 1969) - in 
apparent contradiction to observations. 

In order to understand intra-day variability, several 
intrinsic and extrinsic origins of IDV have been pro- 
posed. The most important extrinsic mechanism is in- 
terstellar scintillation (ISS): turbulent motion of ion- 
ized interstellar matter focuses and de-focuses the light 
from AGN on its way to Earth. ISS has been identified 
observationally by annual modulations of IDV caused 
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by the motion of the Earth around the sun (Gabanyi 
et al. 2007; Liu et al. 2012; Marchili et al. 2012). How- 
ever, the amplitude of interstellar scintillation scales 
with the observation frequency v like v~ 2 - 2 (Rickett, 
Coles & Bourgois 1984). Accordingly, the contribution 
by ISS to AGN variability is substantial at cm-radio 
frequencies but increasingly unimportant at mm-radio 
and shorter wavelengths. In addition, the observed cor- 
relation between rapid variability in radio and optical 
wavelength regimes indicates intrinsic causes (Wagner 
& Witzel 1995). 

The most important intrinsic mechanism for intra- 
day variability is a combination of rapid intrinsic vari- 
ability of the source - especially shocks in AGN jets 
(Blandford & Konigl 1979) - plus Doppler boosting. 
In order to reconcile the observed brightness tempera- 
tures with the inverse Compton limit, Doppler factors 
on the order of 10 to 100 are necessary - which are high 
though still realistic values (Fuhrmann et al. 2008). 

An analysis of plasma-physical processes related to 
AGN outflows usually requires radio observations as 
the emission from AGN jets is - essentially - syn- 
chrotron radiation; the flux density F v approximately 
follows a power law F v oc v~ a with a ~ 0...1. Ad- 
dressing intra-day variability requires flux monitor- 
ing observations that (1) operate at high frequencies 
v > 100 GHz in order to bypass ISS, (2) have good 
(<lh) time resolution, (3) cover observing times as 
long as possible in order to allow for detailed time- 
series analysis, and (4) provide a decent signal-to-noise 
ratio in order to detect even small (<^10%) variations. 
Evidently, these criteria are somewhat in conflict for re- 
alistic observations: on the hand, criteria 2 and 4 sug- 
gest the use of large radio observatories; on the other 
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Table 1. 

Parameter values and corresponding noise levels 



Parameter 


Value 


A 1 Ar. 


u.oo 


At 


100 s 


Av 


1GHz 


T 

± rec 


55 K 


Tatm 


263 K 


^amb 


273 K 


TO 


1.5 


V 


0.9 


PA 


ISOJyK- 1 


T 

J-sys 


290 K 


ST 


0.92 mK 


a N 


0.16Jy 



hand, criterion 3 suggests the exclusive use of a given 
observatory for long periods, potentially many days. 
Accordingly, the scheduling of such programs has been 
notoriously difficult (cf. e.g. Fuhrmann et al. 2008). 

Our work explores an approach complementary to 
observations at large observatories. We analyze the 
performance achievable by using a small, 6-meter class, 
radio observatory dedicated to monitoring of AGN vari- 
ability at 230 GHz. As a realistic example we use the 
Seoul National University Radio Observatory (SRAO) 
located in Seoul, Korea, at the geographic position 
37°27'15"N, 126°57'19"E (Koo et al. 2003; Y.-S. Park, 
priv. comm.). 

2. ANALYSIS 

Our analysis progresses in three steps. First, we esti- 
mate the signal-to-noise ratios (SNR) achievable by an 
SRAO-type observatory. Second, we select potential 
target sources based on AGN properties known from 
the literature. Third, we provide a statistical analysis 
of various types of AGN variability and the correspond- 
ing detection limits. Given the angular resolution of a 
6-m radio telescope, we can safely treat all our targets 
as point sources. Unless stated otherwise, all noise val- 
ues refer to Gaussian la levels. 

2.1 Noise Limits 

The signal-to-noise ratio achieved by a radio observa- 
tory can be calculated like 



SNR = 



Signal 
Noise 



ST 



T^AvAt 



(1) 



ST is the noise temperature of the telescope, Av is the 
observing bandwidth of the receiver, At is the inte- 
gration time, and T sys is the system temperature (cf. 
Wilson et al. 2008). All temperatures are in units of 
Kelvin. For a point source, the equivalent temperature 
Te of the source is simply given by (Thompson et al. 
2004) 



2800 



F v = F v /p A 



(2) 



where k is Boltzmann's constant, and F v is the flux 
density of the source in units of Jansky ( Jy) ; the con- 
version factor pa (in units of Jy K _1 ) is usually referred 
to as antenna efficiency. The effective light collecting 
area A c is given by 



A = A, 



oe 



-(4^C/A) 2 



(3) 



where Aq is the physical aperture area, £ is the root- 
mean-squared difference between actual and ideal an- 
tenna surface, and A is the observing wavelength. For 
the SRAO with an aperture of 6 m, A e /Ao w 0.55 and 
p A = 180 JyK" 1 at 230 GHz. The effective aperture is 
derived from holographic mapping of the antenna sur- 
face (Koo et al. 2003). 

The crucial parameter for our analysis is the system 
temperature T sys which can be written as 



i-sys 



Tree ~t" ^al 



i»7 (l 



- ToXW )+T amh (l- V ) (4) 



where Te is the equivalent temperature of the source, 



where T rcc is the receiver noise temperature, T atm is the 
atmosphere temperature, tq is the optical depth of the 
atmosphere in zenith direction, X{9) is the air mass 
at zenith angle 9, T am b is the ambient temperature, 
and rj is the spillover efficiency including rear spillover, 
scattering, blockage, and ohmic loss efficiency. In the 
present work we use the approximation X{9) w scc(#) 
since we consider only sources with 9 <; 60°. 

In order to arrive at quantitative signal-to-noise ra- 
tios, we make the following choices for the values of the 
various parameters. We assume rj — 0.9, in agreement 
with typical values for modern receivers. For the re- 
ceiver temperature we assume T rcc = 55 K; this is five 
times the quantum noise limit given by hvjk w 11 K 
and typical for modern receivers (cf. Eq. 5.20 of Wil- 
son et al. 2008; h denotes Planck's constant). In agree- 
ment with typical winter conditions at the SRAO site, 
the remaining temperature terms are assumed to be 
Tatm = 263 K and T am b = 273 K, with a - rather con- 
servative - atmospheric optical depth of To = 1.5; to 
take into account an average air mass, we use 9 = 45°. 
Additionally, we assume a bandwidth of Av = 1 GHz 
and select an integration time of At = 100 s. Accord- 
ingly, we find T sys = 290 K and ST = 0.9 mK, cor- 
responding to a noise level in units of flux density of 
ctn = 0.16 Jy. All parameters, as well as the resulting 
noise values, are summarized in Table 1. 
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Table 2. 

Candidate target sources for an SRAO-type radio observatory with SNR>10 



ID (B1950) 


Redshift 


R.A. 


DEC 


F v 


^0 


SNR 


Type 


0059+581°' f 


0.644 


01h02m45.7s 


58d24mlls 


1.9 


230 


12 


LPQ 


0133+476 a 


0.85 


01h36m58.6s 


48d51m29s 


1.7 


94 


11 


HPQ/Blazar/FSRQ 


0235+164 a ' b < f 


0.94 


02h38m38.9s 


16d36m59s 


2.1 


230 


13 


FSRQ/BLLac 


0300+470 


0.47 


03h03m35.2s 


47dl6ml6s 


1.7 


86 


11 


BLLac 


0316+413 


0.017 


03hl9m48.1s 


41d30m42s 


4.8 


230 


32 


Sy2/NLRG 


0355+508 


1.52 


03h59m29.7s 


50d57m05s 


4.2 


230 


27 


LPQ 


0415+379 


0.048 


04hl8m21.3s 


38d01m36s 


4.5 


230 


30 


Syl/BLRG 


0420-014 


0.91 


04h23ml5.8s 


-01d20m33s 


2.5 


231 


16 


Blazar/BLLac/HPQ 


0432+052 


0.03 


04h33mll.ls 


05d21ml6s 


1.6 


230 


10 


LPQ/BLRG 


0528+134 


2.07 


05n30m56.4s 


13d31m55s 


2 


230 


13 


LPQ/FSRQ/Blazar 


0607-157 c 


0.32 


06h09m40.9s 


— 15d42m41s 


2.4 


110 


14 


FSRQ 


0727-115 


1.59 


07h30ml9.1s 


-Ild41ml3s 


2.3 


100 


13 


RLQ 


0748+126 


0.88 


07h50m52.0s 


12d31m05s 


2.3 


230 


14 


LPQ/FSRQ 


0827+243 


0.94 


08h30m52.1s 


24dllm00s 


1.6 


96 


10 


FSRQ/LPQ/Blazar 


0851+202 


0.3 


08h54m48.9s 


20d06m31s 


4 


230 


26 


BLLac 


0923+392 


0.69 


09h27m03.0s 


39d02m21s 


2.6 


230 


17 


Syl/LPQ 


1030+415 


1.11 


10h33m03.7s 


41dl6m06s 


1.5 


94 


10 


HPQ/FSRQ 


1055+018 


0.89 


10h58m29.6s 


01d33m59s 


2.5 


230 


15 


HPQ/BLLac 


1253-055 


0.53 


12h56mll.ls 


-05d47m22s 


9.5 


230 


58 


BLLac/HPQ/FSRQ 


1334-127 


0.53 


13h37m39.8s 


-12h57m25s 


3.8 


230 


22 


HPQ/BLLac 


1510-089 


0.36 


15hl2m50.5s 


-09d06m00s 


1.6 


230 


10 


Syl/HQP 


1606+106 


1.22 


16h08m46.2s 


10d29m08s 


1.8 


94 


11 


FSRQ/LPQ/Blazar 


1633+382 


1.81 


16h35ml5.5s 


38d08m04s 


1.8 


230 


12 


FSRQ/LPQ/Blazar 


1641+399 


0.592 


16h42m58.8s 


39d48m37s 


2.8 


230 


18 


HPQ/FSRQ 


1642+690 a ' b < f 


0.75 


16h42m07.8s 


68d56m40s 


2.3 


230 


14 


HPQ/FSRQ 


1730-130 


0.902 


17h23m02.7s 


-13d04m50s 


2.4 


230 


14 


BLLac/FSRQ/LPQ 


1749+096 


0.322 


17h51m32.8s 


09d39m01s 


2.7 


230 


17 


HPQ/BLLac 


2013+370 




20hl5m28.7s 


37dllm00s 


2.2 


230 


14 


BLLac 


J2044-1043 


0.034 


20h44m09.7s 


-10d43m25s 


2.1 


190 


12 


Syl.2 


2145+067 


0.99 


21h48m05.4s 


06d57m39s 


1.6 


230 


10 


LPQ/FSRQ 


2200+420 d 


0.068 


22h02m43.3s 


42dl6m40s 


1.8 


230 


12 


BLLac 


2216-038 


0.901 


22hl8m52.0s 


-03h35m37s 


1.6 


262 


10 


FSRQ/LPQ 


2223-052 


1.4 


22h25m47.2s 


-04d57m01s 


2.4 


230 


14 


HPQ/BLLac 



NOTES: F v is the source flux at 230GHz in units of Jy, v$ is the observation frequency of data taken from the NED, in 
units of GHz. Entries " — " denote "not available". For some objects IDV has been reported by the studies referred to 
below. 

References: a — Quirrenbach et al. 1992; b — Quirrenbach et al. 2000; c — Kcdziora-Chudczer et al. 2001; d — Kraus 
et al. 2003; e — Lovell et al. 2003; f — Ojha et al. 2004 
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r|=1.000 



50 100 150 200 250 
f 




-0.5- 



-1.0 




■tJa^JkufLu*. jJL_l*^_i— .Ui^^UjiJi ,Lk mAJ~±*A*. ilftkA. 



50 100 150 200 250 
f 





50 100 150 200 250 
f 



Fig. 1. — Illustration of the analysis scheme for periodic variability, a. Idealized, noise-free variable lightcurve with 
amplitude A s = 0.07 Jy and variability frequency / = 5. b. Overlay of 50 Fourier transforms of noise-frec lightcurves with 
A a = 0.07 Jy and / = 5,10, ...,250; as the signals are free of noise, the detection efficiency r/ is unity, c. Like a but with 
Gaussian noise with amplitude ajv = 0.16 Jy. d. Like b but for lightcurves with noise crjv. The horizontal grey line marks 
the 5a detection threshold, e. Like c but with variability amplitude reduced to A s — 0.035 Jy. f. Like d but for lightcurves 
with A s — 0.035 Jy. The substantial decrease in detection efficiency between d and f is obvious. 
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2.2 Source Selection 

In order to select a set of candidate target sources, we 
first extracted a list of 407 radio-bright AGN from the 
NASA/IPAC Extragalactic Database (NED). Unfortu- 
nately, not for all sources data taken at 230 GHz arc 
available. In those cases, we obtained reference fluxes 
F„ , where vq is a frequency as close to 230 GHz as 
possible. We made the usual assumption that the ra- 
dio spectra of our target AGN follow a power law 

F„ = F„ (£f (5) 

where F„ is the source flux at frequency v, F vo is the 
source flux at frequency Vq, and a is the spectral in- 
dex. We assumed a typical spectral index a = 0.6 to 
estimate the source fluxes at 230 GHz. 

We aim at target sources with SNR sufficient to de- 
tect rapid variability. Accordingly, our final sample 
is composed of sources with (1) signal-to- noise ratios 
larger or equal than ten, and (2) transit zenith an- 
gles less than 60°. For each source we computed an 
individual vn using the zenith angle at transit for ref- 
erence and from this the corresponding signal-to-noise 
ratio. For the location of the SRAO, this left us with 33 
sources. The source with the highest SNR turned out to 
be 1253-055 (F 230 ghz « 9.5 Jy, SNR«58), the source 
with the lowest SNR was 1030+415 (F 230 ghz ~ 1.5Jy, 
SNRrsIO). We summarize our sample in Table 2; we 
note that for five of our target AGN intra-day variabil- 
ity has already been reported in the literature. 

We employed various checks to ensure the reliability 
of our sample. Modifying the values for various crit- 
ical parameters, especially a, changes the number of 
selected targets by less than 20%. Therefore we con- 
clude that our sample is fairly insensitive to the choice 
of parameters. 

2.3 Detection Limits of Variability 

Our study also aims at a quantitative description of 
the variability that might be detected under the con- 
ditions outlined above. Accordingly, we analyze sim- 
ulated, though realistic, observations statistically by 
means of Monte Carlo simulations. In the following, 
we always assume a single long - T=14.2h - observa- 
tion of an arbitrary target. The observation time is 
divided into 512 bins of 100 s each. Our calculations 
address two types of potential variability: a periodic 
flux modulation (periodic variability) and a single out- 
burst of activity (flaring variability). Calculations are 
performed in IDL.* 

2.3.1 Periodic Variability 

An artificial lightcurve corresponding to a noisy sig- 
nal with periodic variability and zero average, F s i n , is 



'Interactive Data Language, ITT Exclis Inc., McLean (Virginia) 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 
Variability Flux [Jy] 



Fig. 2. — Detection efficiency r\p of periodic variability 
as function of amplitude of flux variation. This diagram is 
an overplot of 100 analysis runs (grey lines); for each flux 
amplitude, 100 values of rjp are computed, thus showing 
the confidence ranges. 



created by 

F sin = A s sin(27r/i) + F noiso (6) 

where A s is the amplitude of the modulation, / is the 
frequency of the periodic modulation, < t < 1 is the 
time in units of the total observing time T, and F no i se 
is Gaussian random noise with amplitude er/\r. 

We analyze any simulated lightcurve by taking the 
absolute value of its Fourier transform. We classify 
a flux variability as detected if the Fourier trans- 
form has a peak exceeding its average by five times 
its standard deviation.^ In order to derive detec- 
tion thresholds, we examine 50 different amplitudes 
A s = 0.005,0.01, ...,0.25 Jy. For each amplitude, we 
probe 50 different frequencies / = 5, 10, 250. * For a 
given amplitude A s , we define a detection efficiency 

VAA S ) = ^ (7) 

where iV<j p is the number of detections which is 50 
at best (for 50 frequencies of periodic variability). In 



tThis straightforward approach is possible because our artificial 
data are evenly sampled. Realistic flux observations usually re- 
quire the use of periodograms for which the false alarm proba- 
bility follows an exponential distribution (Scarglc 1982). 

t As we normalize the total observing time to the range [0, 1], the 
frequencies arc simply unit-free integers. The Nyquist frequency 
is half the number of data points, i.e. 256 in our case. 
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Fig. 3. — Reduced \ 2 test probing the detection of FLAR- 
ING variability. The reduced \ 2 is shown here as function 
of variability flux Ag, for five different flare time scales oq. 
For each ctq, 100 curves - one for each realization - are 
plotted on top of each other. The horizontal dotted line 
marks the \ /d.o.f. = 1.26 limit above which significant 
variability is detected. 



order to obtain statistically robust results, we test each 
combination of / and A s 100 times, meaning we analyze 
250 000 simulated lightcurvcs in total. We illustrate our 
analysis procedure in Fig. 1 and present the detection 
efficiency as function of amplitude in Fig. 2. 

2.3.2 Flaring Variability 

In order to create an artificial lightcurve with episodic, 
flaring variability, we parametrize the flare as a singu- 
lar flux peak with Gaussian shape. Accordingly, our 
simulated time series is 



Aq exp 



(t-0.5) 2 



2a G 



(8) 



where Aq is the flare amplitude, t is the time in units of 
total observing time, afj is the temporal variance of the 
flux peak (variability time scale) , and F no ; so is Gaussian 
random noise with amplitude <jn- The flare is centered 
at t = 0.5. This parametrization takes into account ex- 
plicitly that the detectability of a flux peak depends on 
its amplitude as well as its time scale. We examine 50 
different flare amplitudes Aq and five flare time scales 
a G = 5/512, 10/512, ...,25/512, meaning 250 parame- 
ter combinations in total. 



Reduced \ 2 Test. Else than for periodic variability, 
there is no obvious recipe for the detection of flaring 
flux modulations. A simple approach is provided by 
a x 2 test that uses a signal constant in time (absence 
of variability) as null hypothesis. Taking into accounts 
the number of degrees of freedom leads to a reduced % 2 
test 



xVcLcf. 



1 



N 
i=l 



F 2 

Gauss, i 



(9) 



where i*bauss,i is the i-th value of the lightcurve de- 
fined by Eq. 8, d.o.f. denotes the number of degrees of 
freedom, and N is the number of data points. In case 
of absence of variability, x 2 /d.o.f . w 1 by construction. 
We calculate 100 realizations of x 2 /d.o.f. for each com- 
bination of Aq and eg, meaning 25 000 calculations in 
total. We examine a range of flare amplitudes Aq = 
0.025,0.05, 1.25 Jy. In Fig. 3 we show x 2 / d -°- f - as 
function of Ac for the five flare time scales examined. 
We consider the null hypothesis "the lightcurve is not 
variable" as rejected if x 2 /d.o.f. > 1.26, correspond- 
ing to a Gaussian significance of 5a for 511 degrees of 
freedom. With increasing flare duration, i.e. increas- 
ing number of data points affected, the minimum flare 
amplitude necessary for detection decreases. 

Multi-Scale Analysis. A more sophisticated but 
nevertheless straightforward approach to detecting flar- 
ing variability is provided by multi-scale analysis (e.g. 
Starck & Murtagh 2006). The underlying scheme is the 
following: 

1. Take a simulated lightcurve and search for a signal 
in excess of five times the noise level of the data. 

2. Bin the lightcurve such that N adjacent data 
points are combined into one averaged value. 

3. Repeat step 1. 

We apply this scheme to our data in an iterative fash- 
ion. In each binning level n = 0,1,2,3,4,5 we com- 
bine two adjacent data points such that J\f=2 n data 
points of the original lightcurve have been combined 
after step n; the case n = simply corresponds to 
the original time series. Binning reduces the noise in 
a lightcurve from on to a^/y/W. Using our usual 5a 
criterion for counting a detection as significant, we are 
able to detect flux peaks with amplitudes 5a n /VN - 
provided the flux peak actually covers approximately 
N data points in the original lightcurve. We illustrate 
our analysis scheme in Fig. 4. 

We analyze 100 realizations for each combination of 
flare amplitude Aq, flare time scale aa, and binning 
level n, meaning 150000 lightcurves in total. We exam- 
ine a range of flare amplitudes Aq = 0.017, 0.034, 0.85 
We define a detection efficiency 



r] F (A G ,a G ,n) = A d , F /100 



(10) 
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0.5 
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(d) 



Binned 
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Fig. 4. — Illustration of the multi-scale analysis scheme for flaring variability, a. Idealized, noise-free lightcurve 
with a Gaussian flux peak with amplitude Ac = 0.35 Jy and width aa = 20/512. b. The lightcurve a with Af=8 adjacent 
data points being binned into one. c. Lightcurve a with Gaussian noise with amplitude on = 0.16 Jy added, d. Binned 
version of lightcurve c. e. Like lightcurve c but with flux peak amplitude reduced to Aq = 0.17 Jy. f. Binned version of 
e. In all diagrams the horizontal dotted lines mark the respective 5a detection thresholds. The impact of binning and flux 
peak amplitude is obvious. 
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Fig. 5. — Detection efficiency tjf for flaring variability as function of flux peak amplitude (variability flux) and binning 
level n derived from multi-scale analysis. Each diagram shows five curves corresponding to the five flare time scales 
indicated in the top left panel. 



where A^j,f is the number of times variability is de- 
tected out of 100 trials for a given set of Aq 1 cg, and n. 
The resulting detection efficiencies are shown in Fig. 5. 

3. DISCUSSION 

Our analysis probes the capabilities of a 6-meter class 
radio telescope with respect to continuum flux obser- 
vations of AGN intra-day variability at high radio fre- 
quencies around 230 GHz. Assuming typical weather 
condition for the location of the SRAO and assuming 
the validity of the usual noise laws of radio astron- 
omy, we conclude that an SRAO-type observatory is 
affected by a statistical (la) noise limit of cat w 0.16 Jy. 
Demanding a signal-to-noise ratio > 10 for AGN flux 
monitoring, we are able to identify 33 sources observ- 
able from the location of the SRAO (Table 2). Thus 
we may conclude that an SRAO-type observatory is in- 
deed able to make valuable contributions to monitoring 
studies of AGN. 

We note that our calculations throughout Sect. 2 
assume somewhat idealized observations leading to 
densely sampled lightcurves without interruptions. In 
realistic observations, the need for occasional calibra- 



tions of flux and amplitude, the observing schemes, and 
technical boundary conditions will result in more com- 
plicated datasets with details depending on the char- 
acteristics of the observatory actually employed. Our 
results depend only weakly on our choice for the to- 
tal observing time T=14.2h: shorter observing times 
reduce the maximum time scale of observable flux vari- 
ations, but not the sensitivity (which is calculated for 
100 s of integration time). 

Unsurprisingly, the limit for the detection of flux 
variations depends on the morphology of the variabil- 
ity. We analyze two rather extreme examples, namely 
periodic and flaring flux modulations. The detection 
efficiency (Figs. 1, 2) of periodic variability - which 
takes into account many different frequencies of mod- 
ulation up to almost the Nyquist limit - approaches 
unity when the amplitude of flux modulations exceeds 
approximately 0.08 Jy. The detection limits for flaring 
variability show a strong interplay between the ampli- 
tude of a flux peak and its duration: shorter flares 
require larger amplitudes in order to be detected. A 
simple x 2 test (Fig. 3) indicates that the minimum am- 
plitudes necessary are approximately located between 
0.3 Jy and 0.8 Jy for characteristic flare time scales 
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gq between 25/512 (in units of total observing time, 
corresponding to «42min) and 5/512 («8min). Us- 
ing a more sophisticated multi-scale analysis and de- 
manding detection efficiencies close to unity (Figs. 4, 
5) reduces these limits to the range from «0.3Jy to 
«0.5Jy. In realistic AGN lightcurves, variability fol- 
lows red noise power laws (e.g. Trippe et al. 2011; Park 
& Trippe 2012) with stronger flux modulations occur- 
ring on longer time scales; red noise variability may 
be regarded as a case intermediate between periodic 
and flaring variability. Notably, we may conclude from 
our analyses that we are able to detect significant vari- 
ability with amplitudes <; 2<jn ~ 0.3 Jy, i.e. on levels 
<J 10% of the average fluxes of our targets. 

The physical conditions within a target AGN can 
be probed via analysis of timescales and amplitudes of 
flux variations. Using realistic astrophysical parame- 
ter values, the observed brightness temperature can be 
expressed (Wagner & Witzel 1995) like 



T b = 4.5 x 10 iu A 



Xd 



t var (l + z) 



K (11) 



where A var is the amplitude of a flux modulation (in 
units of Jy), A is the wavelength (in cm), d is the 
distance to the source (in Mpc), t val is the variabil- 
ity timescale (in days), and z is the redshift of the 
source. For a source at z = 0.3 with A vax = 0.3 Jy 
and t var = 40min, Tb « 10 17 K for observations at 
230 GHz. Even for variability timescales on the or- 
der of one day, T h s=y 10 14 K. Accordingly, any detec- 
tion of intra-day variability with an SRAO-type radio 
observatory corresponds to the observations of sources 
with observed brightness temperatures exceeding the 
Inverse-Compton limit. Using the common assump- 
tion that such an excess in observed brightness tem- 
perature is due to Doppler boosting in AGN outflows 
(e.g. Fuhrmann et al. 2008) one finds an expression for 
the Doppler factor S like 



5=(l + z) 



10 12 K 



l/(3-a) 



(12) 



where z is the redshift of the source and a is the spec- 
tral index (cf. Eq. 5). For observed T b w 10 14 - 17 K, 
5 w 8. ..150 - meaning that we are, a priori, able to 
probe a wide range of plasma-physical and kinematic 
conditions within AGN. 

High-frequency radio monitoring of radio-bright com- 
pact AGN is an important complement to simultaneous 
observations at other wavelengths, both for single-dish 
flux monitoring (e.g. Fuhrmann et al. 2008) as well as 
radio- interferometric mapping (e.g. Lee et al. 2008). 
An important example are measurements of time lags 
between lightcurves obtained at different frequencies, 
providing the dispersion measure - the line-of-sight in- 
tegral of the particle density - and thus information on 
the matter content of the emission region (e.g. Wilson 



et al. 2008). Accordingly, the combination of small, 6- 
m class observatories working at high frequencies with 
larger observatories working at lower frequencies - like 
the Korean VLBI Network (KVN) which operates at 
frequencies of 22, 43, 86, and 129 GHz (Kim et al. 2011; 
Lee et al. 2011) - is likely to provide new insights into 
AGN physics. 

4. CONCLUSIONS 

In this article, we study the feasibility of high-frequency 
radio observations of AGN intra-day variability. We 
assume the use of a small 6-meter class observatory for 
continuum flux monitoring at 230 GHz, using the SNU 
Radio Astronomical Observatory (SRAO) as a realistic 
example. Our work arrives at the following principal 
conclusions: 

1. Assuming an integration time of 100 s per flux 
measurement, a spectral bandwidth of 1 GHz, and 
typical SRAO winter conditions, the statistical la 
noise limit is <tn ~ 0.16 Jy. 

2. Demanding a signal-to-noise ratio > 10, we com- 
pile a list of 33 radio-bright compact AGN with 
fluxes F v ;> 1.5 Jy that are observable at the geo- 
graphic location of the SRAO. 

3. Based on results from exhaustive Monte-Carlo 
simulations, we conclude that it is possible to de- 
tect significant (referring to a 5<r detection thresh- 
old) variability with amplitudes <^ 0.3 Jy, roughly 
10% of the typical fluxes of our candidate tar- 
get sources. Details depend on the morphology 
of variability - periodic, episodic/flaring, or inter- 
mediate. 

4. For nearby AGN (using here z — 0.3 for refer- 
ence), observable flux modulations correspond to 
observed brightness temperatures Tf, w 10 14 - 17 K, 
meaning Doppler factors of 6 w 8. ..150. This in- 
dicates that observations of the type we study are 
able to probe a wide range of plasma-physical and 
kinematic conditions within AGN. 

Overall, we are able to conclude that high-frequency 
radio observations of AGN intra-day variability with 
dedicated small observatories can provide valuable new 
insights into AGN physics. This approach is especially 
powerful when combining those high-frequency radio 
observations with data from other facilities like the 
KVN. 
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